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ABSTRACT 

We consider the pumping of the 63.2/im fine structure line of neutral 01 in the high-redshift in- 
tergalactic medium (IGM), in analogy with the Wouthuysen-Field effect for the 21cm line of cosmic 
HI. We show that the soft UV background at ~ 1300A can affect the population levels, and if a 
significant fraction of the IGM volume is filled with "fossil HII regions" containing neutral OI, then 
this can produce a non-negligible spectral distortion in the cosmic microwave background (CMB). 
01 from redshift z is seen in emission at (1 + z)63.2^m, and between 7 < z < 10 produces a mean 
spectral distortion of the CMB with a y-parameter of y — (10~ 9 — 3 x lCP 8 )(Z/lO -3 Z0)(/j/y ), where 
Z is the mean metallicity of the IGM and Iuv is the UV background at 1300A in units of 1CP 20 
erg/s/Hz/cm 2 /sr. Because O is in charge exchange equilibrium with H, a measurement of this signa- 
ture can trace the metallicity at the end of the dark ages, prior to the completion of cosmic reionization 
and is complementary to cosmological 21cm studies. While future CMB experiments, such as Planck 
could constrain the metallicity to the 10~ 2 Zq level, specifically designed experiments could potentially 
achieve a detection. Fluctuations of the distortion on small angular scale may also be detectable. 
Subject headings: 



1. INTRODUCTION 

The first sources of light that ended the cosmic dark 
ages are expected to have reionized the intergalactic 
medium and polluted i t with metals (see a review by, e.g., 
lLoeb fc Barkanall2001h . To use the metal enrichment as 
a potential probe of the reionization epoch, recent studies 
have concentrated on the neutr al oxygen (OI) pro duced 
by the first massive stars (Oh 2002: iBasu et al.ll2004f ). 
The ionization potential of OI is only 0.2 eV higher than 
that of hydrogen (HI) , and the two species are in charge 
exchange equilibrium. Therefore, oxygen is likely to be 
highly ionized in regions of the IGM where hydrogen is 
ionized. However, the recombination time for oxygen (as 
for hydrogen) is shorter than the Hubble time at z 6, 
indicating that oxygen can be neutral even in regions 
where H has been ionized but where short-lived ioniz- 
ing sourc es have t urned off allowing the region to re- 
combine (|Ohll2002t ). Indeed, \6h~k Haiman I (|2003) show 
that such fossil HII regions can occupy > 50% of the vol- 
ume of the IGM prior to reionization. While the metal 
pollution of these regions is poorly understood, it is fea- 
sible that they contain significant amount of oxygen with 
a neutral fraction close to unity. 

Previous work has proposed to exploit the scattering 
of UV photons by OI, and the corresponding absorp- 
tion features in the spectra of quasars - the OI forest 
(|Ohll2002l ). analogous to the low er-redshift HI Lyman a 
forest (e.g.. iBecker et ain i2006V In the case of HI, an- 
other interesting signature from the high-redshift IGM 
is the 21cm hyperfine structure line, which is made de- 
tectable by UV pumping by the Lyman a background 
(the so-called Wouthuysen-Field effect; see iFieldl 119581 
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or iFurlanetto et al.l 120061 for a recent review) . Here we 
examine whether a similar effect occurs for OI. 

Here we demonstrate that there exist OI lines with the 
required features: the fine structure lines of the electronic 
ground state, 44.1, 63.2 and 145.5/im, can be pumped by 
1300 A photons produced by the first stars or black holes, 
via the Balmer a line of OI. We compute the spectral dis- 
tortion in the CMB created by this effect, and find that it 
can be as high as y ~ 10~ 7 if the OI metallicity at z ~ 7, 
just prior to reionization, is ~ 10~ 25 Z Q . This distor- 
tion could be detectable with future CMB experiments, 
and opens the possibility of performing tomography of 
the reionization epoch using this effect. In combination 
with 21cm studies, it could yield direct measurements of 
the abundance and spatial distribution of metals in the 
high-redshift IGM. 

2. BALMER-a PUMPING OF OI FINE 
STRUCTURE TRANSITIONS 

The basic criteria for a line <-> 1 of a metal species 
or their ions to produce an effect analogous to the 
Wouthuysen-Field pumping of the HI 21cm line, are as 
follows: (i) abundant metals in the IGM in the required 
ionization state; (ii) the <-> 1 transition at a frequency 
suitable for detection, with Einstein coefficients small 
enough to allow the line to deviate from equilibrium with 
the CMB; (iii) the upper and lower states should be con- 
nected via allowed transitions to another state (hereafter 
called state 2), so that they can be "pumped" in a two- 
step process; (iv) a large enough background flux at the 
wavelength A02 ~ A12 corresponding to the <-> 2 and 
1 <-> 2 transitions. In particular, the last criterion im- 
poses the constraint A 2 ^ 1215 A on the wavelength of 
the UV pumping photons before reionization, since neu- 
tral HI depletes the background at shorter wavelengths 
(jHaiman et al.ll2000D . 

We selected oxygen for our study because it can be 
relatively abundant at high redshifts and because the 
similarity of its ionization potential to that of hydro- 
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gen (we leave a systematic search through other metal 
lines to future work). The ground state of neutral oxy- 
gen is split into the three fine structure levels of the 
outermost electrons in the n = 2 shell, hereafter de- 
noted as 0, l a , and l b . We found that among the elec- 
tronic states of oxygen, these are the only one satisfy- 
ing all of the criteria above and which photons would 
be observed in the CMB frequency band. The states 
0-1° and l a -l b are connected by magnetic dipolar tran- 
sitions [01] 3 P 2 -> 3 Pi at 63.2 urn and 3 Pi -> 3 P 
at 145.5 fim respectively and the states 0-l b are con- 
nected by an electric quadrupolar transition 3 P 2 — ► 3 Po 
at 44.1 /im. We shall consider only lines involving the 
state 3 P 2 as the lower level for the transition, since 
the CMB ambient field is practically unable to populate 
the 3 Pi state, and the 145.5 /Ltm line will be suppressed. 
Since these transitions are forbidden, the spontaneous 
emission Einstein coefficients (A^.i^m = 1.34 x 10 -10 , 
A 63 .2 M m = 8.91 x 1CT 5 s-\ Ausj^m = 1-75 x 1CT 5 ) are 
much smaller than the typical values for electric dipole 
transitions (~ 10 s ), a property shared by the hydro- 
gen 21cm line with A H2 \ = 2.85 x 1CT 15 s -1 . 

The states and l a,b are connected to the excited elec- 
tronic state 3 Si in the n = 3 shell (hereafter denoted by 
2) by means of the absorption of an 01 Balmer-a pho- 
ton with wavelength A20 ~ A21 ~ 1300A. A schematic 
level diagram is shown in Figure [TJ In the absence of 
a UV field, 01 is in thermal equilib rium with the CMB 
and, as shown in iBasu et al.l (|2004l ), resonant scattering 
is more important than collisionally-induced emission, 
except at extremely large over-densities. However, if the 
first stars or black holes generated a UV background at 
1300A, then the relative populations of the fine structure 
states and 1 are modified by the two-step "pumping" 
transitions — > 2 — > 1 and 1 — > 2 — > 0. 

For simplicity we start by considering only one fine 
structure transition and only later address the combined 
effect. In what follows we use the indices i,j to denote 
the two levels in the fundamental state with i referring to 
the lowest one. For example for the 63.2.1 /im transition 
corresponds to (0, l a ) (see fig [J). The steady state 
solution for the level population reads 
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This equation defines the spin temperature Ts j%. Here, 
is the equivalent temperature of the i «-> j transi- 
tion (T+ji = huji/ks 



Boltzmann constant), Aji, Bji and Bi 3 



231 (63.2/im/Aji)K with k B the 
are the Einstein 

coefficients for spontaneous emission, induced emission 
and absorption, respectively, PV V and P-j v denote the 
UV de-excitation and excitation rates, l v denotes the 
specific intensity at the resonant frequency of the i <-» j 
transition, and g\a = 3, g^b = 1, go = 5 are the degener- 
acy factors. It is useful to define the UV color tempera- 
ture, Tuvji in terms of the ratio of the UV de-excitation 
and excitation rates, 

P^/PjV = Gfe/Si) expi-T^i/Tuvjt] , (2) 
which represents the spin temperat ure reached in the 
limit of strong pumping. Following [Fill (|1958D . these 
pumping rates can be written as: 
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Fig. 1. — Schematic representation of the transitions of neutral 
oxygen considered here. Note the analogy with the 21 cm transition 
of HI. In the 21 cm transition the relevant line corresponds to a 
hypcrfinc transition, while for OI it corresponds to a fine structure 
transition. The 21cm line is pumped by HI Lyman a photons, 
while the OI line is pumped OI Balmer a photons. 



TABLE 1 

Possible OI transitions between n = 
the n = 3 3 Si level(°'. 
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Wavelength 
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1302.168 


3.41e+08 




1304.858 


2.03e+08 


3 Pl 


1306.029 


6.76c+07 


3 Po 


1641.305 


1.83e+03 




2324.738 


4.61e+00 
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(a) Taken from NIST URL site: 

http: //physics .nist . gov/PhysRef Data/ ASD/lines_f orm. html 
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Here 172 = 3, and the subscript v vq in the right parenthe- 
ses denotes evaluation at v = v pq , i.e. at the frequency of 
the transition connecting states p and q. The sum over m 
takes into account the five possible downward transitions 
from the upper state 2, listed in Table (JTJ. 

From the definition of the UV excitation and de- 
excitation rates, it is easy to see that their ratio is pro- 
portional to the ratio of the quantity c 2 /„/ (2hv 3 ) evalu- 
ated at two slightly different frequencies: v 2 j and i>n. It 
follows then that the color temperature can be approxi- 
mated by 

rp T* 1 T 2l . . 

T uv - Q1 .__ 7H7~i — — = —, > ( 5 ) 



-d\ogn u / dv 



where n v = I v /i> 3 is proportional to the number of 
photons with frequency v, and a s = —d\ogI u /dlogu 
is the logarithmic slope of the background spectrum 
near the OI Balmer a line. We have assumed that 
d\ogn v /dv |„ 2 « d\ogn y /dv |„ 2i . 
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If d log n v jdv |„ 2 . is small and negative then the pump- 
ing mechanism will be more effective in the transition 
i — » 2 than in j — > 2, and as a result the level i will 
be relatively depopulated. In particular, the j — > i pho- 
tons will be seen in emission if Tuv > Tomb, which will 
generally be the case since T 2i ~ 10 5 K 3> Tcmb- 

The amplitude and shape of the UV photon field at 
1300A at the be ginning of rcionization is uncertain. Un- 
like at the HI Ly a frequency of 1215A, however, the 
high-z IGM should be optically thin at 1300A (see dis- 
cussions of the shape o f the soft UV background in 
lHaiman et aTlll996L [2000h . Therefore, the relevant back- 
ground at redshift z should reflect the intrinsic spectrum 
of the dominant sources over the w 27% redshift range 

(1 + z) < (1 + Zsource) < (1 + z) ^ HILyman^/^OIBal 

1.27(1 + z) A . Over the 0.2% fractional difference between 
u 2 j and u 2 i, we then expect a relative decrement of n v of 
~ 0.2%(3 + a s ) throughout vji, or a decrement of order 
~ 1% for sources with a s ~ 2. For our calculations, we 
adopted a fiducial choice of ~ 1%, but in practice, our 
results are insensitive to this value (see below). 

Given the specific intensity /„, the steady state solu- 
tion for the level population provides the following ex- 
pression for the spin temperature Tgji'. 
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z a the presence of dust can 
is given by the CMB, I v 



3. DISTORTIONS ON THE CMB DUE TO 01 

The distortion in the CMB spectrum can be parame- 
terized by the y-parameter, y = AI V / ' B^(Tqmb), where 
AI U is the deviation in surface brightness from the 
Planck spectrum with the unmodified CMB temperature. 
Before computing the distortion using the spin tempera- 
ture above, it is useful to obtain an order-of-magnitudc 
estimate. Because the optical depth of the IGM for the 
UV photons in our case is small (we find t 2 o — 7.3 x 1CU 2 
at redshift 10 for Z/Zq — 10~ 3 ), multiple scatterings will 
be rare. (Note that this is different from the case of the 
HI 21cm, where the Lymana photons responsible for the 
pumping have a large optical depth and scatter multiple 
times.) In this case, roughly a fraction t 2 o of every 1300A 
photon produced by stars will scatter off an 01 atom 
and produce an "excess" vji = vjq photon. Just prior to 
reionization, assuming that a few UV photons have been 
produced per H atom (corresponding to a background 
flux of J v = fewxlO -21 ergs -1 cm~ 2 sr _1 Hz -1 ), the dis- 
tortion should be roughly y ~ fXt^oV ~ 10~ 10 , where 
r\ ~ 6 x 10~ 10 is the baryon-to-photon ratio, and /az 
represents the fraction of CMB photon in the frequency 
band corresponding to the redshift range Az. 

To compute the distortion of the CMB at the redshiftcd 
frequency more precisely, we use the radiative transfer 
equation 

dl v hv 
ds Ait 

where ds denotes the cosmological line element along the 
line of sight. Since the line profile ^{y) is narrow, we 
have ignored the effect of cosmological expansion. After 
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d r„'i\ " t til 1 rv t a ^ r> • ' v taking the 01 abundance in the IGM to be 10 3 times 

& ( £ M *ML t ^™ *** the S0lar value ^adauetalJEOOl, we find the optical 



whereas if A,v <C P 



tion ^ we have T Sji — > i CMB; 
then T SJl -> T uv . 

As long as the slope of the background at 1300A is 
such that as ~ 1, equation ([5]) implies that Tjjv > 
T*ji > 2cmb[^s], so that Tsji must be bracketed be- 
tween Tcmb[.Zs] and Tuv- In the limit of T*j, <C Tuv, 
and thus exp[— T*ji/Tuv] — 1, and Aji 3> Pji (both of 
these conditions are satisfied for the 63.2/xm line), the 
departure of Tg t ji from Tqmb becomes 
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i.e. it depends linearly on the UV flux at 1300 A and 
the dependence on Tjjv is negligible. 

The left panel of Figure [2] shows Ts as a function of 
the UV flux for z s = 10 (solid line) and z s — 7 (dashed 
line) for the case of the 63.2 ^m [0 — » la] line. While Ts 
deviates only one part in 10 3 — 10 from Tomb, this small 
temperature difference nevertheless modifies the popula- 
tion of the levels and 1. 



(9) 

to be very small, (ti<»o = 2.1 x 10 , rii> ~ 3 x 10 13 ). In 
this equation we have approximated the line profile func- 
tion ^(v) by a Dirac delta, and x _1 is the fraction of oxy- 
gen atoms in the 3 P 2 level ground level, (x = 1 to a good 
approximation). Assuming a constant UV background 
(justified by the fact that t 2 q <C 1), the net change to 
the CMB due to the presence of OI can be written as an 
2 i integral along the line of sight of the emissivity, resulting 
in 
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4 Lyman a scattering for these sources should be neglig i- 
ble llChen fc Miralda-Escudl[200l: iFurlanetto fc Pritchard|[2009) . 
while the contribution from yet more distant sources will be 
blocked, because the photons redshifting down to 1300A by red- 
shift z will have passed through the HI Lyman /3 resonance (or 
higher lines), and they will have been broken down into lower en- 
ergy photons. 



(2hvy(c 2 i„) U]0 

The corresponding CMB distortion at the observed 
frequency v — Vjo/(l + z s ) is given by y = 
All /B„(Tcmb[z s ])- In our case, r j0 < 1, and the 
distortion parameter can be written in the the simple 
form y ~ Tj [(expx- l)/(exp (T+jo/Tsjo) - 1) - 1], with 

X = T^.jo/TcMB^s]- 

For the 63.2/im line T^,63.2/jm/7cMB ~ 10 but 
J^,63.2 M m/7cMB(Js, 6 3. 2Aim /7cMB - 1) < 1, thus we con- 
clude that 

,63.2/^m 
Js,63.2fmi \ J- CMB 
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Fig. 2. — (Left) Relative departure of the spin temperature Tg from the CMB temperature as a function of the background intensity at 
the Balmera frequency of OI. Note that these intensities are within the range of amplitudes required for reionization (~ 10 photons per H 

atom, lHaiman eF al, (2000)). (Flight) CMB distortion parameter y = AI^ 3 ' 2 ^" 1 / B v (Tcmb) caused by the emission of OI in the 63.2 /im 
line. In both panels the solid line is for z 3 = 10 an the dashed line is for z s = 7. 



In this case, the distortion parameter y is simply 
the relative deviation of Ts,63.2/xm with respect to 
Tcmb multiplied by the optical depth times the ratio 

2* ) 63.2/im/T's , ,63.2/iin ~ 10. 

The right panel of Figure [2] shows the amplitude of 
the distortion introduced by the OI 63.2 /im transition 
at z s — 7 and 10 as a function of the UV background 
intensity, over a range expected to be relevant for the 
epoch just prior to reionization. Note that, apart from 
the T+fiz ^fim/Tsfis^iim factor, eq. (fTTj) is identical to that 
found for the HI 21cm spin temperature. Since in our 
case, for a metallicity of Z/Zq ~ 10~ 3 the relative de- 
viations of Ts,63.2/im with respect to Tomb are typically 
~ 10 2-3 smaller than in the HI 21cm scenario, and since 
T 2 icm « 10 -4 , our y-distortion is - IO^^Z/IO^Zq)- 1 
smaller than in the 21cm case. For the 44.1 /im (0 — > l b ) 
transition, the amplitudes of y are a factor ~ 5 larger, 
but the relevant frequencies are more contaminated by 
dust and thus more difficult to detect, as we comment 
next. 

4. DISCUSSION AND CONCLUSIONS 

Our basic result above is that the OI 63.2/im (and 
44.1/im) transition could be seen as deviations in the 
present-day CMB black-body spectrum of the order of 
up to 10" 8 at around 160 - 500 GHz (481 GHz corre- 
sponds to z=10 and 160 GHz to z=30, for the 63.2/im). 

The FIRAS experiment has already obtained con- 
straints on the deviation of the CMB spectrum from a 
black-body. We find that in the range of frequencies cor- 
responding t o 10 < z K < 20 f or the 63.2/im line, the 
FIRAS data (|Fixsen et alJ ll996) constrains y < 10~ 5 at 
1-cr level, (4 x 10~ 5 at the 3-o level). Considering the 
expected background UV flux at these redshifts we ob- 
tain a FIRAS constraints of 5 (40) solar metallicity at 
z s > 10. Although this constraint is not yet at an inter- 
esting level, giv en the measurements of metallicities a t 
3 <> z <> 5 (e.g., ISchave et al1l2003t lAguirre et aDl2004T l. 
it is the first direct constraint on the metallicity of the 
IGM in the universe at z ~ 10. The upper limit on the 



metallicity scales linearly with the constraint on y. For 
example, by comparing different, accurately calibrated 
Planck HFI channels constraints y < 5 x 1CP 7 could be 
imposed. An e xperiment able to reduce the FIRAS y un- 
certainty (e.g., iFixsen fe Matherl l2002f ) by two to three 
orders of magnitude could impose interesting constraints. 

The main limitation to the measurement is foreground 
dust emission from the galaxy and from IR galaxies, how- 
ever clean patches of the sky (to remove the dust emission 
from the galaxy) and dust observations at higher frequen- 
cies (to remove the IR galaxies) could be used for this 
purpose. Of course, even if the required S/N is achieved 
in a future observation, there could be other cosmological 
effect s (such as decaying particles; see lFixsen fc Matherl 
2002 for a brief discussion) producing distortion at the 
feeble y ~ 10 -8 level. It remains to be seen, then, 
whether the OI distortion can be disentangled from these. 

One further consideration is that clustering of the 
sources will increase the signal greatly making it more 
easily detectable. The signal should trace the cluster- 
ing properties of the fossil HII regions during reioniza- 
tion and should dominate over the primary CMB angu- 
lar power spectrum in the damping tail. The clustering 
properties of the signal should be qualitatively similar to 
those of the Ostriker-Vishniac effect but will a peculiar 
frequency dependence: the signal will be non-zero only 
for frequencies corresponding to the the 63.2/im transi- 
tion at redshifts with non-negligible OI abundance (and 
UV background). This will be presented elsewhere. 

Our technique complements future experiments that 
will detect the H21cm hyperfine transition in absorption 
because it is sensitive to different systcmatics and op- 
erates at different wavelengths. This would be highly 
valuable, since HI 21cm + OI measurements in combi- 
nation can probe the spatial distribution of metallicity 
directly. 
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